Out-of-equilibrium dynamics of the vortex glass 
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We study the relaxational dynamics of flux lines in high-temperature superconductors with ran- 
dom pinning using Langevin dynamics. At high temperatures the dynamics is stationary and the 
fluctuation dissipation theorem (FDT) holds. At low temperatures the system does not equilibrate 
with its thermal bath: a simple multiplicative aging is found, the FDT is violated and we found that 
an effective temperature characterizes the slow modes of the system. The generic features of the 
evolution - scaling laws - are dictated by the ones of the single elastic line in a random environment. 
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The effect of quenched disorder on the vortex phase 
diagram of superconductors has attracted much inter- 
est since the discovery of high T c cuprates. Since the 
long-range order of the vortex lattice is destroyed by 
any small amount of disorder Q , an amorphous "vortex 
glass" (VG) phase was predicted Q- While it is cur- 
rently understood that there is a "Bragg glass" phase for 
weak disorder there is still some controversy about 
the nature of the amorphous phase at strong disorder. 
Experimentally, one finds an irreversibility line (IRL) be- 
low which the magnetization is irreversible and the linear 
resistivity drops to nearly zero Q . The original proposal 
that the IRL signals a thermodynamic continuous transi- 
tion from the vortex liquid (VL) to a VG with zero linear 
resistivity was initially supported by experiments in 
which a scaling of current-voltage (IV) curves was found 
H and by simulations in randomly frustrated 3D XY 
models with a finite temperature transition 0|. How- 
ever, when magnetic screening is included in these 3D 
XY models the transition disappears Q. Moreover, re- 
cent experiments showed that the IV curves do not scale 
as expected [8J. 

One possible scenario for the IRL is that at low tem- 
peratures experiments occur out of equilibrium. Simula- 
tions of the London-Langevin model for vortices (which 
includes screening) showed that below a crossover tem- 
perature, where the VG-type of criticality is arrested, the 
time scales grow very quickly, similar to a Volger-Fulcher 
form Q as in structural glasses 0. Therefore, it is im- 
portant to compare the dynamic behavior of the vortex 
system with that of structural glasses. In this paper we 
apply the analysis of correlation and response functions 
used in the study of relaxation in structural glasses [Jjj 
to analyze the VG phase in a simulation of a London- 
Langevin model similar to the one studied in Ref. [fillip. 

We consider 3D elastic flux lines in a superconductor 
where the i-th line has coordinates Ti(z) — [xi{z),yi(z)] 
with z the vertical direction along the magnetic field B. 
The superconductor has anisotropy e = £ c /£ab = X a b/X c , 
with £ a (,, £ c the coherence lengths and X a b, X c the pene- 
tration depths; the axis are such that c || z and ab || r^. 



Assuming that r^(z) varies slowly with z, the Hamilto- 
nian for a London model of elastic flux lines is 



Ui(Ar iz ) + U d {r lz ) + ^ U in (r 



i<j 



Having discretized z in units of d z , = r^(z) rep- 
resents the 2D position of a line element in the plane 
z. The interaction energy between line elements in the 
plane z is approximated as J7j n = 2eod z K (r / X a b) with 
e = ($ /47rA a6 ) 2 0,0 03. The elastic line energy of 
the i-th vortex is Ui = \ci(d z Yi) 2 dz, which is discretized 
as Ui = \ci{^) 2 d z , with Ar lz = r %tZ+ i - r i>z . A natu- 
ral choice for d z is the distance between CuO planes, in 
which case one can account for the effect of the Josephson 
coupling using XJ\ = 2c;Aj-^^ for |Ar| > 2Xj (keep- 
ing the previous expression for |Ar| < 2Aj), where Xj = 
d z /e is the Josephson length. The quenched disorder po- 
tential due to impurities is t/<j(r) = J d 2 r'u(r')p(\r — r'|), 
where p(r) = 2£ 2 fc /(r 2 + 2£ 2 6 ) and (ii(r,#V)) = 
jS(r — r')8 zz > defines the disorder strength 7 We 
model the dynamics with the Langevin equation 



dv lz {t) SH[{r lz (t)}} 



V- 



dt 



Sri 



where 77 is the Bardeen-Stephen friction coefficient. 
The thermal force if z {t) satisfies (/£,u(t)) = 0) an d 
</£,J*)/^>(*')> = 2 v k B TS(t - t')5^5 U 'S^, where 
fj,, fi , = x,y and T is the thermal bath temperature. 

The above model gives a good quantitative descrip- 
tion of the vortex physics of moderately anisotropic high- 
T c superconductors like YBa 2 Cu 3 07_ 5 |H H Q [H . 
We therefore choose parameter values corresponding to 
YBa 2 Cu 3 7 -5: e = 1/5, X ab /U = 100, and XjJU = 
16; and we use a = e 2 e 2 [1 + \n{\ ab /d z )] /it U\. The 
strength of disorder is set to 7 = 10 -5 , for which case we 
find that above B cr ~ 0.002iJ C 2 the Bragg peaks dissa- 
pear and the flux lines are frozen in a highly amorphous 
structure at low temperatures. We therefore choose to 
study the case with B = 0.01iJ C 2 ^> B cr , which is deep 
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FIG. 1: (a) Dynamic wandering and (b) mean square dis- 
placement for elastic lines with no in-plane interaction and 
without disorder, after the quench to T — 0.06. The elastic 
correlation times t\ and th separate three regimes (single pan- 
cake, elastic line and center of mass diffusion) . The behaviors 
B ~ (t - t w ) and B ~ (t - t w ) 1/2 are highlighted. 

within the VG regime at low T. Time is normalized 
by to = Cab^Aoi length by the vortex lattice parameter 
a = [2<i>o/(v / 3S)] 1 / 2 , energy by €od z , and temperature 
by eodz/ks- We simulate N — 56 vortices in a box of 
size 7ao x 8aoV3/2 with periodic boundary conditions for 
the in-plane coordinates. The z direction is discretized in 
L = 50 planes with free boundary conditions. Averages 
are performed over 10 realizations of the disorder. 

In order to study the out-of-equilibrium dynamics we 
use the following protocol. First we equilibrate the sys- 
tem at an initial high temperature well inside the VL 
(Ti = 0.3) evolving during t — 10 4 steps. Then we quench 
the system to a low temperature T, where the time count 
is set to zero. Starting with this far from equilibrium ini- 
tial condition, the system is let evolve during a waiting 
time t w , after which the quantities of interest are mea- 
sured. In general, one defines two-time correlation or 
response functions C(t,t w ). When the system reaches 
equilibrium, these quantities become independent of t w 
and depend only on the difference r = t — t w . If the 
system is not able to reach equilibrium within the obser- 
vation time window, C(t,t w ) depends on the two times. 
In particular, if the decay gets slower for longer t w 's we 
say that the system "ages" . We study: (a) The dynamic 
wandering jig, W(z, t, t w ), defined by 

NW(z,t,t w ) = ^<|[r lz (<)-r l0 (i)]-[r J2 (^)-r i0 (^)]| 2 ) 

i 

which measures how the displacement of the vortex seg- 
ment in the z-th plane with respect to the bottom plane 
correlates between t and t w ((. . .) means average over 



thermal noise and disorder), (b) The mean square dis- 
placement (MSD) in the planes, 

B(t,t w ) = -^X)<[*<*(t) -XiS w )f). 

iz 

Firstly, we analyze the single, non-interacting, flux line 
without disorder. In Fig. ^a) and (b) we show the dy- 
namic wandering and the MSD, respectively, when the 
line is quenched to T = 0.06. In this case the flux 
line reaches equilibrium in a short time: both correla- 
tion functions are independent of t w and depend only on 
t = t — t w (data for several t w overlap in the plots). 
From the behavior of W, wc define a set of charac- 
teristic times t z , such that when r = t — t w > t Zl 
W(z 7 t » t z ) ~ z 2t > saturates to a constant value. Here Q 
is the roughness exponent given by thermal fluctuations, 
i.e. C = Ct = 1/2 .143 - t z is the time needed for the 
line element in the 2-th plane to feel the elastic interac- 
tion with the line element in the plane z = 0. In Fig. ^ 
we show the times scales t\ for neighboring planes sepa- 
ration and tL for full system size separation. We found 
that for r < ti, W(z,t) ~ r for all z. For times t > t^, 
W(L,t) saturates. For times t\ < t < t z the dynamic 
wandering shows an intermediate regime between these 
two extreme behaviors. In Fig. db) the MSD is shown 
and the same regimes are identified. First, for r < t\ the 
MSD follows the 2D diffusion of individual line elements 
("pancakes"), B(t) ~ r, before the elastic interplane in- 
teraction becomes relevant. Second, for t\ < r < ti,, 
flux line thermal relaxation is observed, characterized by 
sublinear diffusion B(t) - r a , with a = 1/2 [H llif. 
Third, for r > t^, we observe diffusion of the center of 
mass of the flux line B(t) ~ r. is the time scale above 
which finite size effects dominate. The time regime we 
wish to study is the one without finite size effects, t < 
(t < 10 4 for L = 50). 

Secondly, we analyze the evolution of interacting lines 
in the presence of disorder. When we quench the sys- 
tem to relatively high temperatures (T > 0.2) it reaches 
equilibrium, W(t, t w ) and B(t,t w ) are independent of t w 
and they behave as in Fig. ^ for non-interacting flux 
lines meaning that we are clearly within the VL. Be- 
low a crossover temperature T g w 0.18 the system is 
no longer able to equilibrate, and the correlation func- 
tions depend on t w . The MSD at T — 0.02 is shown in 
Fig-IU a )- For t < t\, the dynamics is governed by single 
pancake fluctuations and the MSD is independent of t w . 
For t > ti, there is "aging", the longer the waiting time 
t w , the slower the relaxation of the flux lines 0] . 

In order to study the modifications of the fluctua- 
tion dissipation theorem (FDT) in the out-of-equilibrium 
regime, a response function should be measured. To this 
end, a random force of the form fj Z — (5si Z x is switched 
on at a time t w on a replica of the system, where 5 is the 
intensity of the perturbation, and Si Z = ±1 with equal 
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probability 0, . The integrated response is 



where xf z and Xi Z correspond to the position evaluated 
in two replicas of the system, with and without the per- 
turbation. In equilibrium, FDT implies 2T\(t) = B(r). 
In Fig.H», we show 2T X (t,t w ) at T — 0.02. X depends 
on t w showing aging. It is also clear that for long r and 
long t w 2Tx is not proportional to B, violating FDT. 
This type of behavior: aging in B and x an d violation of 
FDT for long t w and r is observed at all T < T g . 

Typically, a correlation function C(t,t w ) in structural 
glasses has an additive scaling form C(t,t w ) = C eq (t — 
t w ) + C„.„\h(t) / h(t ,„)] , with, very often, h(t) = t (known 
as simple aging) [ill l2(i| . This does not hold in the 
vortex problem. Instead, we found a "multiplicative" 
scaling similar to the one proposed for the low temper- 
ature behavior of a directed polymer in random media 
|2lj | and the out-of-equilibrium critical dynamics of the 
2d XY model H3- Following Yoshino j2l| we tried the 
scaling form B(t,t w ) = B{t)t* and 2T X (t,t w ) = x(t)t w , 
with t = t/t w and B and x given by 



B(t) 

m 



Cl (T)(f-l) Q ( T ) f«l 

c 2 (T)(f-l) Q ( T ) t»l 

Cl (T)(f-l) Q ( T ) f«l 

y(T)c 2 (T)(t - i»l 



where Cj, c 2 and y are temperature dependent coeffi- 
cients |2J. y{T) measures the modification of FDT I20J 
2T X {t,'t w ) = y(T)m,t w ) (or x = y(T)B), and an effec- 
tive temperature |2J] is defined by T e g =T/y. 



FIG. 3: (a) Multiplicative scaling for the mean square dis- 
placement and the integrated response function correspond- 
ing to the data in Fig. [5] T — 0.02, and the value of a is 
quoted, (b) The parametric plot x(-B) for the same data. 
FDT holds at short rescaled times while a violation of FDT 
with y = 0.111 holds at longer times. The inset shows the 
short rescaled time FDT regime for t w = 10 4 . 



In Fig. EI a) we show the scaled B and £ at T = 0.02. 
The data for different t w fall on two master curves. B 
and x coincide for t <C 1, which means that FDT holds, 
while for longer times (? > 1) y(T) < 1 signals the vi- 
olation of FDT. In Fig. Ofb) a parametric plot of the 
scaled B and x is shown, for t <C 1 it has a slope equal 
to one, corresponding to the FDT regime (see an enlarg- 
ment in the inset), while for longer times the slope is 
y = 0.111. From this kind of plot the temperature de- 
pendent dynamic exponent a(T) and the parameter y(T) 
in Fig. are obtained. At high temperatures a = 1/2 
as for single flux lines. The low temperature phase has a 
lower dynamic exponent (cv(T) < 0.4), implying a much 
slower relaxation. In Fig. 4(a) we see that a(T) depends 
weakly with T and that it decreases with increasing dis- 
order strength within the glassy regime. In Fig. E^b) 
we show y(T) = T/T e g; it is well described by a lin- 
ear form implying an effective temperature T e g that is 
independent of T. A linear fit yields T c s — 0.175. It 
is remarkable that this value is very near the crossover 
temperature T g « 0.18 below which aging is observed. A 
similar result is observed in structural glasses: T c ff as T„ 



(as in a random energy model scenario 18J ). 

To further check a "fragile glass transition scenario" 
|l0| . we estimated a characteristic relaxation time by 
evaluating C(t,t w ) — exp[— B(t, t w )]. In the VL we ex- 
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pie violation of FDT with an effective temperature T e g 
which is independent of T within our numerical accuracy. 
No, because aging follows a multiplicative scaling, sim- 
ilar to polymers in random media and, more precisely, 
critical systems like the 2d XY model 22] . Longer times 
are needed to decide whether this "critical" scaling holds 
in the limit t w — ► oo and to grasp the nature of the low- 
T phase. Our results suggest to do experiments with a 
fast quench to low T: with relaxational measurements of 
transport or magnetic properties one could study the ag- 
ing regime while with voltage or flux noise measurements 
one could test how the FDT is violated. 

We thank discussions with H. Yoshino and T. Grig- 
era and support from SECYT-ECOS, Conicet, CNEA, 
ANPCYT and ICTP grant NET-61, and Fundacion An- 
torchas (SB). LFC is a member of IUF. 
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FIG. 4: (a) Temperature dependence of the dynamic expo- 
nent. The equilibrium value for single elastic lines, a = 1/2, 
is highlighted, (b) The slope y(T). The data is fitted linearly 
with y — T/T a ff and the values of T B g are indicated. Inset 
shows the relaxation time at the VL-VG transition. 



pect C(r) ~ exp[— (r/t r ) a \. For single flux lines the In- 
dependence of t r should be t r oc 1/T 2 0,0]. We ana- 
lyzed this correlation function for different T, fitting the 
corresponding stretched exponential with a = 1/2 and 
estimating the relaxation time t r (T), that is plotted in 
the inset of Fig.^Ja). In the VL we obtained t r oc 1/T 2 , 
as expected. Below T g « 0.18 the system does not equi- 
librate but we can extract a by fitting C(t,t w ) for 
the largest t w with the same exponential form. The ob- 
tained t r increases very rapidly below T g , reflecting that 
the system falls out of equilibrium in this aging regime. 

Finally, we performed the same analysis for single flux 
lines in the presence of disorder. For high temperatures 
we found no aging and a w 1/2, with important finite 
size effects [cfr. Fig. [J. Below a crossover temperature 
we observe similar aging to the one shown in Fig. [3 The 
multiplicative scaling, and the same definitions of B and 
X, describe the data very accurately, the only differences 
being that a(T) takes larger values that are closer to 
a = 1/2 than in the VG [see Fig. Ufa)], and that T cS 
takes a smaller value, T e g = 0.156, than in the VG. 

The results described in the last paragraph imply that, 
rather surprisingly, the out-of-equilibrium slow dynamics 
of the VG is dominated by the relaxation of the elastic 
lines along the z-direction. Aging follows a multiplicative 
scaling as for single lines, but with a smaller exponent 
a(T). Do vortices freeze like a structural (fragile) glass? 
Yes and No. Yes, because we found aging and a sim- 



[1] A. I. Larkin and Yu. N. Ovchinnikov, Zh. Eksp. Teor. 

Fiz. 65, 1704 (1973) [Sov. Phys. JETP 38, 854 (1974)]. 
[2] D. S. Fisher et ai, Phys. Rev. B 43, 130 (1991). 
[3] T. Giamarchi, P. Le Doussal, Phys. Rev. B 55, 6577 

(1997). 

[4] K. A. Muller et al, Phys. Rev. Lett. 58, 1143 (1987); Y. 

Yeshurun and A. P. Malozemoff, ibid. 60, 2202 (1988). 
[5] R. H. Koch et al, Phys. Rev. Lett. 63, 1511 (1989); A. 

M. Petrean et al, Phys. Rev. Lett. 84, 5852 (2000). 
[6] J. D. Reger et al, Phys. Rev. B 44, R7147 (1991); T. 

Olson and A. P. Young, Phys. Rev. B 61, 12467 (2000). 
[7] H. S. Bokil and A. P. Young, Phys. Rev. Lett. 74, 3021 

(1995). 

[8] D. R. Strachan et al, Phys. Rev. Lett. 87, 067007 (2001); 
I. L. Landau and H. R. Ott, Phys. Rev. B 65, 064511 
(2002). 

[9] C. Reichhardt et al. Phys. Rev. Lett. 84, 1994 (2000). 
[10] M. D. Ediger et al. J. Chem. Phys. 100, 14200 (2000). 
[11] L. F. Cugliandolo in Slow Relaxations and Nonequilib- 

rium Dynamics in Condensed Matter, J. -L. Barrat et al. 

Eds., (Springer, Berlin, 2002). 
[12] A. van Otterlo et ai, Phys. Rev. Lett. 81, 1497 (1998). 
[13] S. Ryu and D. Stroud, Phys. Rev. B 54, 1320 (1996). 
[14] G. Blatter et al, Rev. Mod. Phys. 66, 1973 (1994). 
[15] In stronly anisotropic systems like Bi2Sr2CaCu2 0s+<s a 

model of interacting 2D pancakes [J. R. Clem, Phys. Rev. 

B 43, 7837 (1991)] may be more adequate. 
[16] H. Yoshino (unpublished). 

[17] Aging was also found in a 2d model by M. Nicodemi and 
H. J. Jensen, Phys. Rev. B 65, 144517 (2002). 

[18] G. Parisi, Phys. Rev. Lett. 79, 3660 (1997). 

[19] A. B. Kolton et al, Phys. Rev. Lett. 89, 227001 (2002). 

[20] For a review see A. Crisanti and F. Ritort, J. Phys. A: 
Math. Gen. 36, 181 (2003). 

[21] H. Yoshino, Phys. Rev. Lett. 81, 1493 (1998). 

[22] L. Berthier et al, J. Phys. A: Math. Gen. 34, 1805 (2000). 

[23] We assumed a linear dependence of x with B, but a non 
trivial function x(B), as found in the spin- wave approx- 
imation to the 2d XY model |2^|. cannot be discarded. 

[24] L. F. Cugliandolo et al, Phys. Rev. E 55, 3898 (1997). 



